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Abstract Reaction-diffusion equation model for a 2- 
dimensional gas discharge system is introduced in relationship 
with the recent experiment by Nasuno 00 [2]. The model shows 
formation of spots, molecule-like organization of a cluster of 
spots, moving of spots, and intermittent collapse and recre- 
ation of spots, in agreement with experiments. The pattern 
dynamics are classified into distributed, localized, and moving 
spots, and periodic wave pattern. The phase diagram displayed 
against the current and pressure has some agreement with that 
observed in the experiment. Dynamic change of spot numbers 
as well as the structure of the spot cluster is studied in terms 
of the collective variable of charge, and discussed as chaotic 
itinerancy [3]. 

1 Introduction 

Pattern dynamics in non-equilibrium systems have 
been studied over decades, in fluid, chemical reaction- 
diffusion systems, liquid crystal, and so forth. The 
extensive studies in the field have elucidated a rich 
variety of pattern dynamics, together with the ad- 
vances in theoretical analysis [I] ~ [T5] , 

One field in the pattern dynamics that is not so 
well explored in comparison with the above exam- 
ples, is a discharge system. When a strong voltage is 
applied between the electrodes, electric discharge ap- 
pears through the gas filled in the chamber. The dis- 
charge often forms a complex spatiotemporal pattern, 
as has been studied theoretically and experimentally 

(CSMZzl). 

In such discharge system with a variable resister, 
there exists some global constraint among each lo- 
cal discharge processes. Viewed as pattern dynamics, 
this means existence of global coupling among local 
dynamical processes. On the other hand, global cou- 
pling among nonlinear units often shows a non-trivial 
collective motion, as has been extensively studied as 
collective motion in globally coupled dynamical sys- 
tems. Now in the discharge system, there exists in- 
terplay between collective motion in globally coupled 
systems and pattern dynamics by local nonlinear dy- 
namics. As a study of nonlinear system, it is inter- 
esting to search for some novel dynamic state, due to 



this interplay. 

Indeed, there is a beautiful experiment by Nasuno, 
suggesting such novel, non-trivial dynamics [T] [2J. He 
set up an experiment consisting of two parallel plates 
between which electric charge flows. By controlling 
the voltage and current, he found formation of spots, 
organization of 'molecular- like' structures of spots, 
complex motions of them, including the dynamics 
what he called teleportation. 

In the present paper we introduce a phenomenolog- 
ical model on 2-dimensional (2D) gas discharge sys- 
tem, as a coupled electric circuit subjecting a global 
constraint, to describe the salient features observed 
in Nasuno's experiment, and to make further predic- 
tions. As a theoretical model, we adopt a two-layer 
model that describes well a glow discharge system 
[19] ~ [24]. This two-layer model consists of a non- 
linear resistive part and a linear resistive part. By 
extending this two-layer model, we construct a cou- 
pled dynamical system consisting of elementary cir- 
cuit systems, each of which represents the process 
at the interior between square electrodes. By taking 
spatial continuum limit, this model is reduced to a 
two-dimensional reaction diffusion (RD) system with 
global constraint. Through extensive simulations of 
the model, we will show several phases of the pat- 
tern dynamics, and reproduce some salient features 
in Nasuno's experiment. 

On the one hand, our model is an abstract and 
simplified model, and it may not completely corre- 
spond to the discharge system by Nasuno. On the 
other hand, the model shows several novel interesting 
patten dynamics as a reaction-diffusion system with 
global coupling, which itself deserves investigation. 

The outline of the paper is as follows. In Sec [2] 
the experiment by Nasuno is briefly described. In 
Sec [3] we present some assumptions in order to re- 
produce the experiment and then construct a cou- 
pled circuit system as a model on the gas discharge 
experiment. In Sec [4] results of extensive simulations 
are presented, with classification of distinctive pat- 
tern dynamics as in the phase diagram. Formation of 
spots, and a cluster of spots are numerically found, 
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as well as the dynamics of creation and annihilation 
of spots, while the dynamical systems mechanism of 
the process is discussed. In Sec|5] summary, discus- 
sion and future perspectives are presented. 
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Figure 1 : Schematic figure of 2 layers model for glow 
discharge. We consider that the glow discharge struc- 
ture consists of two parts, the non-linear resistive re- 
gion and the linear one. 



2 Experiment 

Here we briefly describe a remarkable experiment by 
Nasuno In the experiment, two square elec- 

trodes are immersed in low pressure air P and con- 
nected to a dc power source e and a variable resistor 
R ex . The power source e can supply either a constant 
voltage of up to 2 kV or a constant current of up to 
150 mA via a series resistor Rex 35.2 kf2. 

He adopted the parameter region of the gas dis- 
charge experiment so that the product of pressure 
P and distance d between electrodes satisfies P x 
d ^0.04 (Torr-cm), by using Air. The parameter 
region is located in the vicinity of the minimum of 
paschen curve, according to Fig. 7. 3 of chapter7 in 
[28) . It is known that, under the experimental condi- 
tion, current-voltage characteristics has monotonous 
curve j25]~[3D], and that no pattern dynamics is usu- 
ally observed. He fixed the total current I instead of 



voltage, which is different from the conventional dis- 
charge experiments. However, when the current I is 
fixed as a constant, V is determined uniquely by a 
monotonous current-voltage characteristic so that V 
is also constant. Indeed, according to the experiment, 
current- voltage characteristic I— V is monotonous, so 
that I and V remain fixed [2 . However, local current 
and charge density between electrodes can change in 
space and time, to produce nontrivial dynamics. 

By controlling I and P, complex pattern dynamics 
in glow discharge region was observed. By increasing 
current I, the pattern dynamics of discharge changes 
as follows. 

• Just after onset of glow discharge, one isolated 
spot with a typical size (that is, a humped bell- 
like light intensity distribution) appears. The 
spot wanders irregularly on the square plane. 

• With the increase of 7, the number of isolated 
and wandering spots increases. 

• As I is increased further, some spots form a 
molecule-like localized structure, which is called 
cluster. This cluster also wanders through the 
plane. 

• With the further increase of J, the number of 
spots in the cluster increases leading to various 
clusters with a variety of configurations of spots. 
Switching among different clusters occurs inter- 
mittently. 

• At much higher I, the excited domain (the area 
with high light intensity) forms a string or a 
closed loop. 

• A remarkable phenomenon which he termed 
'Teleportation' was reported, where annihilation 
of a spot and immediate recreation of a different 
spot at a distant place was observed. For exam- 
ple, when 3 spots existed, one of them suddenly 
was disappeared, and then a new spot immedi- 
ately appeared at a distant position, located in 
the neighbor of a remaining spot, so that the 3 
spot-state is recovered PQ. 

The phase diagram for the above behaviors is dis- 
played in terms of I and P in the paper [1] . 



3 Model 

3.1 Assumption 

Here we introduce an electric circuit model corre- 
sponding to the above experiment. The model is 
essentially based on a circuit model introduced by 
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Purwins et al.[21]~[24], derived from an analogy be- 
tween an electric circuit and a reduced equation from 
plasma physics. The model is at a macroscopic 
phenomenological level, where discharge process be- 
tween 2D electrodes is described as charge transfer 
process by a coupled circuit system . For each ele- 
ment, we adopt the so-called "2 layer model" which 
consists of non-linear resistive part and linear re- 
sistive one, as schematically shown in Figjl] The 
model has succeeded in explaining several discharge 
experiments |19j ~ |24j . Here we revise their model so 
that it meets with the experimental condition by Na- 
suno. 

To set up the coupled circuit model, the following 
assumptions are made following Purwins et al. 

• The experimental condition is set at the region 
of glow discharge. 

• We describe the phenomena in terms of current 
density i(x, y) and space charge density q(x,y) 
because we focus only on the charge transfer pro- 
cess between two electrodes. 

• Gas pressure effect between electrodes known ex- 
perimentally is taken into account. According 
to [23] [21], this effect is effectively described as 
diffusion of current density that flows between 
electrodes. Based on the experimental results, 
we assume a monotonous relation between gas 
pressure and the diffusion of current density. 

In considering Nasuno's experiment, we make fur- 
ther assumptions: 

• Global current-voltage characteristic shows a 
monotonous cubic curve in a logarithmic scale 
[31] . as will be given later(Fig{3](b)). 

• The total current that flows into electrodes is 
kept as constant . Therefore, the voltage drop 
between electrodes also remains constant. 

• In electrodes (represented by 2-dimensional 
planes), local voltage-current density character- 
istic is assumed to be cubic in a logarithmic scale 
with a negative slope. 

3.2 Specific Model 

The first set of assumptions underlie the coupled cir- 
cuit model by Purwins et al. [5TJ ~ [53] , while we ar- 
range it to meet with the latter set of assumptions to 
adapt to Nasuno's experiment. 

The discharge model which we consider here con- 
sists of two square electrodes (area S) and the inte- 
rior, which are under "External circuit" that consists 
of the dc power source e and the variable resistor 



R ex , as shown in Fig(2] The total current I flow- 
ing into the electrodes is controlled to be constant fi a . 
Thus V is uniquely determined to be a constant. Fur- 
ther, each elementary circuit k is connected with the 
neighboring ones by linear resistor R and is also con- 
tacted with "External circuit". As shown in Fig[2] 
this elementary circuit consists of a linear resistor, 
capacitance, linear coil and a nonlinear resistor given 
by a specific i-v characteristic. The linear resistor r, 
the region U and the intermediate region correspond 
to linear resistive layer, the non-linear resistive layer, 
and the interface in 2-layer model, respectively. Each 
part is characterized by a linear resistor r, capaci- 
tance C, the constant self-inductance I of the coil, 
while the voltage drop v k at the non-linear layer fol- 
lows the i-v characteristic given by a log-scaled cubic 
function of i}.. The total current I is distributed to 
each clement with a current Jk, while its current at 
B and U is given by jk and ik + q\ respectively, as 
shown in Fig(2] 

Now, the main discharge system is described by a 
set of TV ordinary differential equations as N coupled 
circuits, as shown in Fig [2] In the model it is assumed 
that there are some charge flow into the interface (i.e. 
at "B" in the figure) from the non-linear resistive 
layer. This flow of charge density is represented by s'k. 
The voltage drop at the interface is represented by 
Since the total space charge Q should be constant in 
time between electrodes, there is a global constraint 



on J2k s k as w iU be discussed again in Sec 3.3 



Considering the law of the conservation of charge 
(the first Kirchoff 's rule) at the point A, we get 

Jk - ik + + s k) - (qk-i + Sfc-i) (1) 

+(qk + Sk) - (qk+i + Sfc+i)}] = 0, 

while the law of the conservation of charge at the 
point B gives 



dq k , ds k 
]k - {lk+ — ) = — 



(2) 



At the region U, the voltage balance (the second Kir- 
choff 's rule) leads to the following equation 



,di k . qk „ 

l lt +Vk -c =0 - 



(3) 



In addition to (jTJ ([3j) , we need to consider the con- 
straint on the constant total current / and constant 
voltage drop V. 



N 

I = Jk = const. 

fc=i 

V = rJ k + H + = const. 



(4) 
(5) 



Here, to solve the above equations, we assume a 
constraint for the description of Sk as a function of 
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{i ra , q m }, (m — 1,2, N). Accordingly we must de- 
fine Sk satisfying Q^Q. After some calculations, 
with Vk = v(ik) and q' k = cjk + Sfc, Sfc are determined 
so as to satisfy 

dt r 



1'k 


ik H 


1 


rC 




di k _ 




Sk 


dt ~ 


IC 


IC 



N 



I = ik = const. 

k=l 

N I 

Q' = Y.<lk = Cv(-)N = const. 



fc=l 



•N' 

i Q! 



V = r- + -^ = const. 



We define Sk as 



Sk 



N 



fc=i 



(Sec appendix for details.) By taking spatial contin- 
uum limit, we obtain the following 2-component RD 
equation, with D q > = -^q,Q' = / q'dS, as 



dq 1 _ V 
~dt ~ ~ 
di 



^-i + D q ,Aq' 
r rC 

4 s v(i) 

di = ic~ic~~r + 1 

s= l -{Q' -C fv(i)dS} 



I 



idS = const. 



Q' = j JdS = Cv (^) S = const 
V = r— + v(-) = const. 
By using dimensionless variables 



(6a) 
(6b) 
(6c) 
(6d) 
(6e) 
(6f) 



where r= ^,a = r^,£>'= ^ q ,D= ^, 
£ 2 = rCD q (characteristic length), S = |y. 
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Figure 2: Schematic figure of our 2-dimensional dis- 
charge model. Our model system consists of N cou- 
pled circuits. The element k is connected with neigh- 
bor ones by linear resistor R and with "External cir- 
cuit" . Although displayed in a one-dimensional rep- 
resentation, what we actually simulated in the paper 
is a two-dimensional case. 



g 

Civ 



I. V 



v, q^q,t 



the RD equation (6a I 



di 
dt 



dq , r 

T— = V 

dt 

= ~{q-s- 
a 



I 



{Q 



t = -f and redefining 

|6fj ) is written as 
- q ~ ai + Aq 
v(i)} + DAi 
v(i)dS} 



V 



I = J idS = const. 

qdS = v( — )S — const. 
T 

const., 



a 7s + v{ s ) 



3.3 Flowing charge density 



(7a) 
(7b) 



Here we make some remarks on the flowing charge Sk 
at the point B in Fig [2} The point B corresponds to 
the interface of the 2-layer model, that is, the posi- 
tion between the "Negative glow" and Faraday dark 
space, which exists for the steady glow discharge (Fig 
(7c) [Tj. i n this region, electrons flowing from Cathode 
layer (non-linear resistive layer) combine with posi- 
(7d) tive ions flowing from "Positive column" , or anode 
(i.e., the linear resistive layer). There, some complex 
(7e) processes occur through diffusion, excitation, ioniza- 
tion, and recombination of molecules. Hence there is 
a longitudinal flow of charge, and the flowing charge 
Sk at B is distributed. 



(7f) 
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For a homogeneous steady glow discharge, there is 
no charge in this region. Furthermore, Sk at B is 
zero, as the initial state before the discharge. Hence, 
if a spatially homogeneous state were stable, the flow 
term Sk would remain to be 0. However, as will be 
shown, for almost all the parameter regions, such ho- 
mogeneous state is unstable, where charge density q 
is distributed inhomogeneously in the non-linear re- 
sistive layer. In this case, flowing charge density s is 
also distributed in the interface. Hence we need to 
consider this term. 

After taking a spatial continuum limit, the total 
flowing charge is given by the spatial integration, 



h= / sdS(=Q- / v(i)dS), 



(8) 



where h means the total charge flowing from the 
non-linear resistive layer into the interface, which is 
a macroscopic dynamical variable characterizing the 
global charge transfer process between the non-linear 
layer and the interface. 

In Nasuno's experiment [T] , [5] , due to the short dis- 
tance between electrodes, the Positive column does 
not appear. Hence the anode and the anode surface 
correspond to the linear part and the interface, re- 
spectively. Since the anode has very small resistance 
in the experiment, s is assumed to be distributed over 
the anode surface. 

3.4 Parameter setting 

Hereafter we study the behaviors of the present model 
given by the equations (7al~d7f||, by controlling the 



parameters / and D, which are the total current into 
the system, and a parameter for the gas pressure ef- 
fect, respectively. Hence the control parameters I 
and D correspond to those in the experiment, i.e., 
the total current and the pressure P, respectively. 

As initial conditions of i(x,y) and q(x,y), we 
mostly choose "uniform state i u , q u " with small fluc- 
tuations r](x, y) as a uniform random number over 
[-0.01,0.01]. Accordingly total current / (= i u S), to- 
tal charge Q(= J qdS) and voltage V are determined 
initially (See appendix). The phenomena to be dis- 
cussed is always observed from these initial condi- 
tions, i.e. the pattern dynamics is attracted to a 
global attractor. Unless otherwise mentioned, the 
other parameters are fixed as 

r = 20, £ 2 = 0.004 (characteristic length), a = 
400^ = 0.0108,S"= A = 250. 

Numerical simulations are carried out by the mesh 
size 128 x 128. Unless otherwise mentioned, we choose 
a periodic boundary condition. 



These parameters are chosen so that a single el- 
ement satisfies the behavior for a simple discharge 
(without spatial pattern). Before describing the pat- 
tern dynamics of the model, we first show how a single 
element behaves. Fig [3] is the nullcline of the single 
element dynamics (Figj3Ta)), and global I—V charac- 
teristic (Figj3Tb)). We mainly focus on the case that 
the steady state (i u ,q u ) is among the critical point 
(1,1) in Figga) m- 




Figure 3: Nullcline of the single element dynamics 
(a), and global I—V characteristic (b). We mainly fo- 
cus on the case that the steady state (i u , q u ) is among 
the critical point (1,1), where it is stable if i u < 1, 
whereas destabilized if i u > 1 and then the element 
shows limit cycle oscillation. I—V characteristic is 
monotonous so that a steady state is realized. Q is 
also a function of i u . v(i) is cubic and log-scaled 
function of i. Note the relation / = i u S. 



4 Phenomena 

4.1 Global phase diagram 

According to our interpretation of the present model 
([7a|^([7f|, the discharge occurs at local site if 
i{x,y) > i u . Hence we show the pattern dynamics 
by displaying only the pixels that satisfy i{x, y) > i u , 
an active region with discharge, by gray. 

Here we study the pattern dynamics by changing 
iu{= j;) and D, by focusing on the parameter region 
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with i u ~ 1, where the discharge phenomena occur. 
As a 2-dimensional phase diagram with regard to i u 
and D, the pattern dynamics we have observed is 
classified into five phases, as displayed in Fig|4j We 
first give a brief description of each phase, and will 
discuss the characteristic of each phase in detail later. 

For small i u (^ 1.0), there are three phases, de- 
pending on the value of D, that are "distributed spot 
phase (DS)", "localized spots phase (LS)", and mov- 
ing spot phase (MS), respectively, as shown in Fig[4j 
In these phases, spots are formed. Discharge occurs 
locally, only within these spots. Indeed, when i(x, y) 
is plotted in space and time, it shows a stepwise in- 
crease at the border of each spot, and shows a high 
plateau within each spot. The pattern dynamics here 
will be discussed from the configuration of spot pat- 
terns and their motion, with which the three phases 
are classified. For all these three phases, the number 
of spots increases with i u . For a larger value of i u , 
there appears a periodic pattern in space, which we 
call Periodic Wave (PW) phase. For a further large 
value of i u , there is a uniform glow (UG), without 
spatial inhomogcneity. The behavior of each phase is 
summarized as follows; 

• DS phase — Spots are arranged with some dis- 
tance on the average, which decreases with the 
increase of i u . For small i u , a few spots are iso- 
lated with some distance, while hexagonal pat- 
tern of spots appears as i u is increased. After 
transient time, these spots are fixed, and do not 
move. 

• LS phase — Spots form several clusters similar 
to molecular structures, as reported in the exper- 
iment by Nasuno. The clusters show a dynamic, 
complex behavior with the formation, collapse, 
and regeneration. 

• MS phase — A spot can move by itself in the 
space, which shows a soliton-like behavior with- 
out collapse. 

• PW phase — There appears a periodic wave pat- 
tern in space and time. 

• UG phase There is a uniform glow, without 
spatial inhomogeneity. 



Collective variable — J sdS (= h) — As a global 
measure characterizing the pattern dynamics, the in- 
tegration of s, h = / sdS is often useful. In the 
subsequent subsections, the pattern dynamics will be 
characterized by the motion of h, which is a collective 
variable. Now we discuss the pattern dynamics, by 
referring to the change of this collective variable h. 
(Note that h -C Q in the present parameter setting.) 
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Figure 4: Phase diagram as function of i u (= ^) and 
D. Vertical axis is log-scaled. Global property of 
the model is roughly classified into five phases under 
i u > 0.9. 



4.2 DS phase 

Isolated steady spots Starting with initial con- 
dition i(x,y) — i u + rj(x,y), competition for current 
among elements occurs, resulting in instability of the 
homogeneous state. This leads to the formation of 
spots with localized high current. These spots are of 
the same size. Initially, each spot moves very slowly, 
and is arranged so that each is located with an equal 
distance, until a regular, stationary spot pattern is 
formed. 



Hexagonal pattern As i u is increased, the num- 
ber of spots is increased, so that the distance between 
spots is decreased. As i u is increased from 1.5 to 2.0, 
a regular hexagonal pattern of spots is formed, as 
shown in Fig(5] The formation of this regular lat- 
tice of spots is understood through weak repulsion 
between neighboring spots. 

Behaviors at the boundary between DS and 
other phases As i u is increased toward the bound- 
ary to the PW phase, the spots start to breathe, and 
their sizes (and shape) change, while the distance be- 
tween spots still remains almost equal as in Fig{6]ja- 

1) . The collective variable h changes intermittently, 
corresponding to irregular pattern dynamics (Fig{6]ja- 

2) )- 

Around the boundary between the DS and LS 
phases, the characteristic wavelength disappears, and 
the configuration of spots is irregular, while spots 
start to show a bursting behavior, as well as split- 
ting, as in (Fig |6jb-l)). The change of h as well as 
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Figure 5: Steady patterns observed in the DS phase. 
Spots are of the same size under fixed D. Initially, 
each spot moves very slowly, and is arranged so that it 
is located with the equal distance. As i u is increased, 
the number of spots is increased. As the number is 
increased with the increase of i u , the distance be- 
tween spots is decreased. Left-Isolated steady spots. 
i u =0.87,0.93,1.0 in the order from the top to bot- 
tom with, D=0.01. Right-Hexagonal spots pattern. 
i u =1.2,1.5,2.0 in the order from the top to bottom 
with D=0.01. The mesh size here is 256x256 
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Figure 6: Behaviors at the boundary between the 
DS and PW phases(a) and at the boundary between 
the DS and LS phases (b). The left figures (1) give 
snapshots of pattern dynamics, while the time se- 
ries of the integration of flowing charge h(— J sdS) 
corresponding to them are plotted in (2) (right). In 
(a), the spots start to breathe, so that their sizes 
change, while the distance still remains almost equal. 
In (b) , the distance is no longer equal, and an irregu- 
lar pattern is formed. (a) i u = 2.4, D = 0.01, (b) 
i u = 2.0, D = 0.012, mesh size 256x256 



the pattern dynamics is irregular(Fig[6jb-2)). 

In the experiment by Nasuno, the behavior corre- 
sponding to this phase is not observed. We expect 
that the range of the change of pressure in the ex- 
periment may not be sufficient to cover low values 
needed and to detect this phase. 



4.3 LS phase 

For small i„(~ 0.9), a single spot appears. As i u is 
increased, the number of spots increases, through suc- 
cessive split of the original spot. Spots are not com- 
pletely stable in time, and some of them collapse after 
breathing. As replication and extinction of spots are 
repeated, a cluster of spots is formed. For small i u , 
the number of spots changes between 1 and 2, while 
the range of the number of spots that the system can 
take increases with i u . Some of the generation pro- 
cesses are displayed in Fig [9] 



Phase diagram within the LS phase By mea- 
suring the dependence of the maximal number of 
spots on i u and D, the phase diagram within the 
LS phase is depicted as shown in FigjTJ where the 
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maximal number of spots is plotted in the D — i u pa- 
rameter space. Hereafter denote the case "Mk " when 
k spots exist at maximum. On the other hand, if n 
spots exist at a moment, we call the state "sp n ". In- 
deed, the configuration of the phases in the diagram 
of FigjT] obtained from our model agree rather well 
with that obtained experimentally by NasunopQ. 
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Figure 7: Phase diagram of the LS phase. The be- 
haviors in the LS phase are classified in terms of the 
maximal number of spots, and plotted as a phase di- 
agram in the D — i u parameter space. Note that, in 
this phase diagram, "M n " shows the case in which 
the maximal number of spots is n. 



Localized cluster of spots; molecule-like struc- 
ture In the following, we discuss dynamics of a lo- 
calized cluster of spots. After spots are generated by 
splitting, they separate up to some distance. Thus 
they form a chain-like structure, a localized struc- 
ture like a molecule. As the distance between spots 
is increased and a spot is separated from a cluster 
very slowly, it starts breathing. Following the ampli- 
fication of this oscillation, the spot collapses. After 
this collapse, split of some other spot(s) follows, and 
the number of spots returns to the maximal under 
the given value of i u . In this way, a cluster of spots 
is sustained. 

When the maximal number of spots > 6, there ap- 
pears a variety of configurations for the cluster, as 
shown in Fig [8] Depending on which spot splits to 
which direction, there appears a different configura- 
tion. Over a long time span, one can observe a variety 
of spot configurations. 

In a long time scale, a cluster wanders slowly on a 
2-dimensional space, since the configuration after re- 
construction of a cluster is slightly different from the 
original, even if they are similar. After the split of a 
spot, there remains some asymmetry between the two 
spots, which causes a wandering motion of a cluster. 
We show three examples of such wandering cluster, 
by displaying successive snapshot patterns in Fig. WU 
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Figure 8: Snapshots of typical clusters, sp2 ~ sp$. 
They are the snapshots of the clusters with the max- 
imal number of spots under the given i u and D. 
£>=0.014. sp 2 ««=0.93,sp 3 i„=0.934,s]54(a)i u =0.98, 
sp 4 (b)i u =1.01, sp 5 (a)£ u =l.l,sp s (b)i u =1.06, sp 6 (a) 
i u =l.l, sp 6 (b)i u =l.l,sp 6 (c)i u =l.l, sp e (d)i u =l.l, 
sp 8 (a)i u =1.15,sp 6 (b)i„=1.15. 
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Figure 9: Generation processes of the clusters sp3 
sps through successive split of the spots, SP3 : spi - 



sp 2 sp 2 ■ 
0.014). sp 4 
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> sp3 => sp3 cluster (i u — 0.939, D — 
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D = 0.014). sp 5 : sp 2 
sp 5 cluster (i u = 1.06, D = 0.013). 
sp 2 — ► sp 4 sj»4 — > sp 6 =>■ sp 6 cluster 
1.08,1? = 0.013). sp 7 : sp 2 - 
spi spi cluster (i u = 1.1, D 
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spi 

0.013). 
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1.1, D 
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sps 



sps cluster 



where the maximum number of spots is 2 (M2) for 
Fig[l0];a), 3 (Af 3 ) for Fig^b) and 4 (M 4 ) for Fig,[l0 
(c), respectively. The motion of cluster is not unidi- 
rectional, but its direction changes irregularly. This 
irregular motion arises since the two spots after split 
are not completely identical. 



Dynamics of the collective variable h corre- 
sponding to the change in the spot number 

When the number of spots changes, the transfer of 
charge is altered drastically. Hence, it is relevant 
to study the number change in relationship with the 
collective variable h, i.e., the integration of the flow- 
ing charge. In Fig [TT|a),(b), and (c), we show the 
time series h for M2, M3 and M4 cases. We plot 
the time series of h (in (a-l),(b-l),(c-l)), as well 
as the orbit in the phase space, by embedding the 
time series of h into three dimensional phase space 
(h(t),h(t + l),h(t + 2)). 

For all the cases, h increases with the spot number 
n. When the spot number stays at some value, h 
remains almost constant. Note that h value is mainly 
determined by the spot number of the moment sp nj 
and is almost independent of the maximal number 
Mk for a given condition. For example, h for sp2 
under M2 almost equals to that of sp2 under M4. If 
the system stays at a state with given sp n (n < N), 
it is clearly seen as each plateau in the time series 
of h, or a region with residence of an orbit around a 
fixed point in the phase space. With the breathing 
of spots, h also starts to oscillate, and in the phase 
space picture, the orbit spirals out of a fixed point, 
leading to a switch to a state with a different number 
of spots. 

Consider the case M3 or M4. There the decrease in 
the number of spots occurs as (sp4 —>)sp3 — * sp2 — * 
sp\, successively. During this decrease, the system 
stays at each state with an intermediate spot num- 
ber, for some time interval, while the return process 
spi — > sp2 — > sp 3 (or SP4), is relatively rapid, without 
staying long at each intermediate state. (The state 
sp\ is unstable, and the orbit exits immediately). 
These decrease and increase in the spot number are 
repeated. We also note that the switching process ob- 
served here is true for a state with a higher number 
of spots, as have been numerically confirmed up to 
sp s . In Fig[l3j we display an example for M 6 ~ M s 
case. 

In the three-dimensional representation of the 
phase space from the time series of h, each state of 
a given spot number seems to be regarded as a sad- 
dle which has one-dimensional stable manifold, and 
two-dimensional unstable manifold, with an unsta- 
ble focus. The attraction and repulsion of each orbit 
around a saddle appears whenever the orbit passes 
through a given sp n state. Each state of a given spot 
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number is approached from a certain direction, and 
the orbit spirals out from it. This process is analo- 
gous to that observed in Shilnikov chaos In the 
Shilnikov chaos, however, the spiral motion is stable, 
giving a stable 2-dimensional manifold for a focus, 
while a third direction gives an unstable manifold. 
In contrast, in the present case, the former is unsta- 
ble, and the third direction gives a stable manifold. It 
is expected that the instability here leads to irregular 
wandering of spots. 

However, when h increases, the orbits path through 
several saddle points with three real eigenvalues in 
this representation of the phase space, as long as sp n 
is less than its maximal spot number (as can be seen 
in Fig. 



11) 



Note that the orbits approaching the 
saddle and leaving it take different paths, depend- 
ing on either if h is increased or decreased, as shown 
in Fig 12 In other words, the three-dimensional 



phase-space representation is insufficient, and the 
high-dimensionality in the original dynamics leads 
to the dependence on the history of the change in 
h. In this sense, the itinerant motion over differ- 
ent spot numbers will be better described as chaotic 
itinerancy [34J, [35], in the sense that several low di- 
mensional ordered states are visited through high- 
dimensional chaotic motion. 

Transition rules between several clusters Here 
we study transition-rules between clusters in some 
detail. 

As the increase in the spot number occurs through 
the splitting of a spot, there exists a certain transi- 
tion rule from a state with a lower number of spots to 
that with a higher number. Here we display the paths 
of generating sp n (n — 2, ...,8) through the splitting 
process (see Fig|T4|. As shown, each cluster with 
a larger number of spots is systematically generated 
from a specific state with a smaller number. Note 
that there are several states with different configura- 
tion of spots for a given number of them, when it is 
larger than or equal to 6 |37j . 

This transition rule of clusters is enriched with the 
increase of i u . In Fig[l5] Schematic figure for the 
change of the transition rule with i u is displayed, 
corresponding to the pattern in Fig 14 where the 
solid upward and broken downward arrows show mul- 
tiplication of spots by splitting and their collapse by 
breathing, respectively, while the thickness of the ar- 
rows shows the frequency of such processes observed. 
The transition rules are summarized as follows. 
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Figure 10: A sequence of snapshots for sp2, sp3, sp^ 
state. We see that the positions of sp2, sps, sp4 shift 
in time due to the asymmetry for splitting. In a 
long time scale, the clusters appear to wander, (a) 
sP2 : i«=0.93(M 2 case), (b) sp3:i„=0.939(M 3 case), (c) 
sp 4 :i«=Q.98(M 4 case), L>=0.014 



• When i u is small, only the transition sp\ 
occurs. 



SP2 



• As i u is increased, the states sp3 and sp^ (which 
is one type of the 4-spot state) appear, with the 
transitions s Pl ~ sp 3 s Pl ~ spf 
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With the further increase of i u , another 4-spot 
cluster with asymmetric configuration appears, 
Now there are transitions 



denoted by sp^'. 

spi ~ sp 3 , spi 



(1) 

sp\ , 



spi 



(2) 

sp\ ■ 



As i u is increased further, the states sp^ ~ sp$ 
appear, with a variety of transitions. These tran- 
sitions are divided into two groups as in Fig[i4| 
that is, those among spi ~ spi 1 and among 

(2) 

spi ~ sp\'. 



Note that with the increase in the spot number, 
both the configurations and transitions are diversi- 
5 fied, while h shows complex dynamics corresponding 
''* 1 \o the diversification (Fig 13 1. 
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Figure 11: Dynamics of h. (1) Its time 
(2) orbits embedded into a 3-dimensional 
The case Mi (i„=0.93), (b) The case M 3 ( 
(c) The case M4 (i u =0.98). Each plateau 
corresponds to each state spi ~ sp^ with 
different i u under fixed D. D = 0.014 



When the maximum number of spots is larger 
than or equal to 4, there exist clusters with the 
same spot number but different configurations. 
Those states give almost the same value of h, 
and they are degenerated in the representation 
of ft. 

When the maximum number of spots is larger 
than or equal to 6, the collapse process is di- 
versified, and indeed is more diverse than that 



series and 
phase, (a) 
i„=0.939), 
value of h 
respect to 



displayed in Fig 15 
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Figure 12: Motion around a saddle in high dimen- 
sional phase space. The nature of the saddle point 
sp m is different depending on if h increases or de- 
creases. The approach to and deviation from a 
saddle changes according to the in(de)crease of h. 
i u = 0.939, D = 0.014 (M 3 case, m=2) 



The higher the symmetry in the configuration 
of the clusters is, the higher is the frequency of 
the appearance of such configuration. Indeed, 
this frequent appearance of configuration with 
higher symmetry is also observed in the experi- 
ment ([50). 

The transitions between the degenerate states, 
i.e., different configuration states with the same 
number of spots, cannot occur directly. Only 
after the increase or decrease in the spot number, 
they can mutually change, as in Fig. |15| 



Correspondence with the experiment The 

pattern dynamics we observed here agree with that 
observed experimentally [T] [5] as is summarized as 
follows. 

• The global phase diagram on the spot number 
with regards to i u and D (pressure) agree rather 
well. 

• The cluster of spots forms a molecule-like struc- 
ture, whose structures are identical. These struc- 
tures are similar to that observed on the experi- 
ment. (In the experiment, the molecular state is 
so far reported up to the spot number 6 mainly) . 
In our case, the number of spots shows intermit- 
tent decrease, while by comparing the state that 
is dominant in time, the agreement is clear. 
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Figure 13: Time series of h for i u =l.l, 15=0.013. h 
shows complex dynamics, corresponding to the di- 
versification of the cluster configurations with the in- 
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Figure 15: Schematic figure for transitions among 
quasi-stationary states. The solid upward and bro- 
ken downward arrows show multiplication processes 
of spots by splitting and extinction process of them by 
breathing, respectively. The thickness of the arrows 
indicates the frequency of the processes. The typical 
cluster configurations in the present figure are shown 
in Figflll 



• When the maximal spot number is larger than or 
equal to 6, a variety of forms of clusters appear, 
both in our model and in experiment. 

• An interpretation of teleportation in the term 
of Nasuno: the reported phenomena in which 
one of the two spots disappears, and right af- 
ter it, a new spot appears around the remaining 
spot. Since it is observed as if one of the spots 
moved to a distant place within a very short time 
scale, he called this phenomenon as 'teleporta- 
tion'. According to our result, we can give a 
possible interpretation to this phenomenon. In 
the M2 regime, one of the spots disappears inter- 
mittently as we have mentioned. Following this 
disappearance of one spot, the remaining spot 
divides into two. Note that the time interval be- 
tween the disappearance of one spot and the di- 
vision of the remaining spot is very short. Then, 
within the resolution of experimental measure- 
ment, this process is observed as if one spot tele- 
portated to the location close to another spot. 



Figure 14: Paths to create state with a higher number 
of spots, sp n (n=2,...,8) by splitting. Note that only 
the states generated by splitting are displayed, and 
some other complicated ones, which are generated as 
a result of extinction of spots by breathing, are not 
displayed here. 



4.4 Moving spot phase 

At this phase a single spot moves in the space. For 
0.93 < i u < 1.0 in Fig[4] a single moving spot appears, 
while multiple moving spots are observed with the 
increase of i u . 
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Figure 16: Dependence of the orbits of a single mov- 
ing spot on the value of i u , under periodic boundary 
condition where a solid line shows the center of mass 
of the spot. We show the direction of movement of 
msl with an arrow. D =0.027, mesh size 256x256. 



Soliton-like spot Here a single spot (abbreviated 
by "msx " ) moves without split or collapse. The lo- 



cus of a spot is displayed in Fig [T6j with a periodic 
boundary condition, where a solid line shows the cen- 
ter of mass of ms\. The center of mass is defined by 
the mean position among sites with i(x,y) > i u . Of 
course, if the spot is completely symmetric in shape, 
it cannot move. The motion of a spot is due to asym- 
metry in the shape of ms\, and indeed the spot is 
deviated from a circle slightly. 

The nature of motion changes with the increase of 
i u . For small % u value, that is just above the onset 
of the appearance of msi, the motion is linear with 
a constant speed. (See Fig 16 i u = 0.925). With the 



increase of i u , the motion starts to have a curvature, 
and the locus is bended. This curvature is increased 
with i u , and the locus shows a circle, as in Fig \TE[ 
i u = 0.934, 0.94. With the increase of curvature, the 
speed also increases. 

For this motion of spot, the choice of boundary con- 
dition may be more crucial. For example, to compare 
with an experiment, the periodic boundary condition 
may not be relevant. Hence, we have also made some 
simulations with Neumann boundary condition, as 
shown in Fig [17] With reflection at the boundary, 
the motion becomes more complex. 



iu=D. 945 






Figure 17: The loci of a single moving spot under the 
Neumann boundary condition. Arrows show the time 
course of the spot center. Because the spot reflects 
at the boundary, the motion is complex as shown for 
i u = 0.94. On the other hand, it often shows periodic 
orbits as given for the loci for i u — 0.93,0.935,0.945. 
For i u = 0.95, 0.97, it shows almost circular motion, 
but the center of the circular orbits shifts gradually. 
These behaviors depend on the ratio of the radius of 
curvature of the spot locus to the electrode size . 
D =0.027 



Multiple moving spots As i u is increased further, 
a single moving spot splits, leading to moving mul- 
tiple spots(Fig 18 1. Here, the moving spot is first 
distorted to extend its tail as spiral. When the tail is 
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3) 



Figure 18: Successive snapshots of multiple moving 
spots. We show the direction of movement of the 
moving spots with arrows. (1)A single moving spot 
just before splitting. Here the size starts to expand. 
(2) It splits into two spots, but one of them is too 
small and vanishes. The other distorted spot sur- 
vives and soon begins splitting again. (3) Just before 
splitting of the distorted spot. (4) Just after splitting 
of the distorted spot, leading to the state ms-i. Each 
spot is separated from each other in the opposite di- 
rection. (5)The state ms2 is maintained, but one of 
them starts breathing and will vanish before long. 
i u = 1.05,.D = 0.03 



Jt 



Figure 19: Banana-like spot. A spot is spread in the 
direction perpendicular to the direction of its mo- 
tion, and finally the banana-like shape is formed. 
The banana-like spot moves straight without col- 
lapse. i u = 1.2, D — 0.03 



Coll 
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Figure 20: (a) Periodic pattern. The parameter val- 
ues are i u = 5.0, D = 0.015, while the similar pat- 
terns appear over a wide range of parameters, (b) 
Plane wave. i u = 5.0, D = 0.03. This plane wave 
moves in one direction as shown by the arrow. 



small, the spot rotates with it, while as it is larger, 
the spot is divided into two, which moves to the oppo- 
site direction (generation 0/TOS2). After the motion 
of these two spots, collapse of one spot occurs as in 
the LS phase, and the system comes back to a single 
spot state. These processes are repeated to lead to a 
complex motion of spots. 

As shown in Fig 18 the spot shape starts to be 
distorted. This distortion is much clearer, as in is 
further increased, where the system approaches the 
boundary with the PW phase. There the spot size 
is larger, and forms a distorted " banana- like" shape 
(Fig 19 1. This " banana- like" spot moves straight 
without collapse. 

Periodic wave phase In general, as i u is in- 
creased, each site tends to synchronize with each 
other, since the effective (attractive) coupling among 
oscillator elements is stronger. Accordingly, in this 
regime of a large current, spots no longer exist, and 
they are replaced by a pattern periodic both in space 
and time (FigpOja)). 



Traveling wave One typical pattern observed here 
is plane wave(Fig |20|b)). As shown, an extended 
string of discharged region propagates with a con- 
stant speed here. A string of concentrated electric 
discharge is formed over the plate, and it travels with 
a constant speed. 



Stripe pattern Here we briefly describe pattern 
observed at a different set of parameter values of r 
and £ than that adopted so far. For some parameter 
values of t and £, we have found a steady pattern, 
which is different from the DS phase. As an example, 
we discuss briefly the pattern observed for r = 35 
and £=0.008. As i u is increased, the connected spots 
form a string and form a stripe pattern as shown in 
Fig [2T] The length of string is increased with i u , 
and finally a stripe pattern covers the whole space 
as in Fig|2T] Through the time evolution, a steady 
pattern is formed, without temporal change. Such 
pattern formation is also common to that observed 
in Bcnard convection with a large aspect ratio. 
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Figure 21: Steady stripe patten. Spots are extended 
and connected to form a string (i u — 2.2,4.2). As 
i u /*, many strings are formed, and as a result they 
are connected with the neighboring strings as shown 
(stripe). Plotted for i u = 6.2,8.2. 



5 Summary and Discussion 
5.1 Summary 

To sum up we have introduced a coupled circuit 



model (7a I ~ (7/1, corresponding to the discharge 



experiment by Nasuno. Mathematically, the model 
belongs to a class of reaction-diffusion system with 
global coupling due to the global constraint on the 
conservation of charge and current. When the to- 
tal current is large, the system shows a homogeneous 
glow, while this homogeneous state is destabilized 
with the decrease in the current, and a variety of pat- 
tern dynamics of discharge is observed, as are classi- 
fied into four phases, depending on the value of total 
current and the diffusion constant that corresponds 
to the pressure value in the discharge experiment. 

For the regime with a lower current, the regions 
with concentrated discharge form spots. They are 
classified as 

• Distributed Spot phase: Spots are arranged 
with some distance, to form a regular array of 
spots. The phase appears for small D. 

• Localized Spot phase: A few number of local 
spots exist, which form a cluster. With the in- 
crease in the current the spot number increases. 
In the cluster spots are arranged like a molecule, 
while there are several configurations of spots 
when their number is larger than or equal to 4. 
Also, we have observed intermittent collapse of 
spots to decrease their number, and immediate 
recovery by the splitting of spots. These pro- 
cesses are repeated as a cycle, and with this cy- 



cle, the cluster of spots wanders throughout the 
space. 

• Moving Spot phase: 

For a larger value of I?, a shape of a spot is asym- 
metric, and each spot starts to move by itself. 
With the increase of the total current the mo- 
tion changes from linear to circular, while com- 
plex motion is observed when the boundary con- 
dition is not periodic. 

Besides these spot phases, we have observed an- 
other phase at a high current region: 

• Periodic Wave phase, where a string of dis- 
charged region propagates in space. 

These pattern dynamics are characterized by intro- 
ducing the collective variable h — J sdS, the inte- 
gration of the flowing charge s/., which expresses the 
global charge transfer between the non-linear resis- 
tive region and the interface in the direction of the 
gap. Indeed, h turns out to characterize the number 
of spots, and the birth-and-death dynamics of spots 
are represented by the temporal change of h. 

5.2 Comparison with the experiment 
by Nasuno 

The behaviors observed in the LS phase reproduce 
the phenomena observed in experiments. They in- 
clude: 

• Increase of spot number with the current. 

• Molecule-like structures of spots and their vari- 
ety in shapes. 

• Disappearance and recreation of spots, as is 
termed 'Teleportation' by Nasuno. 

• Wandering motion of spots. 

loop In the experiment, at a high current region, 
a loop structure that moves in space is observed. Al- 
though a connected string observed in the PW phase 
is similar to the loops observed in the experiment, 
one difference here is that in our case the strings are 
connected through the whole space. Since we have 
chosen the periodic boundary condition, the edges of 
the traveling string in Fig [20] are connected, to form 
a loop, but such traveling string observed is extended 
to the whole space in our model. 

So far, we are not yet confident if the strings in the 
PW phase correspond to the loops in the experiment. 
As the loops are not small in contrast to localized 
spot structures, the influence of the boundary may be 
important. A suitable choice of boundary condition 
should be necessary. Instead of periodic or Neumann 
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boundary condition, the Dirichlet boundary may be 
more relevant to make more accurate correspondence 
with the experiment. This problem is left for the 
future. 

Related studies 

There are some recent studies concerning Nasuno's 
experiment on gas discharge. 

Static pattern of particle structure is discussed by 
modifying Gray-Scott equation in [39], while its re- 
lationship with discharge system is not so straight- 
foward. In |38| . pattern dynamics caused by local 
current heating of the gas are studied by using the 
left-branch of Paschen curve, termed as obstructed 
discharge. In our study, we consider surface charge 
on the electrode in the obstructed discharge in Sub- 
section [373] and study instability of the homogeneous 
steady state at the left-branch of Paschen curve. 

5.3 Discussion 

Novel phenomena in reaction-diffusion sys- 
tem Since our model belongs to a class of reaction- 
diffusion systems, many of the pattern dynamics ob- 
served here are common with those studied therein; 
formation of spots, array of spots, and strings. Still, 
the molecule structure of spot cluster, cyclic pro- 
cess of collapse and split of spots, wandering mo- 
tion of spots are rather novel and characteristic to 
the present model. For these phenomena, inclusion 
of global coupling into local reaction-diffusion equa- 
tions is essential. Systems both with local and global 
couplings have also been studied in surface catalytic 
reactions [40], coupled maps [41]. and so forth. Search 
for a novel class of pattern dynamics by global cou- 
pling will be interesting in the future. 

Itinerancy over clusters with different config- 
urations of spots Through split and collapse of 
spots, their number changes, and also transitions be- 
tween quasi-stationary states with different clusters 
take place. This transition among quasi-stationary 
states is found to be governed by a specific rule with 
regard to the increase or decrease in the spot num- 
ber. With the aid of the phase-space representation 
by the macroscopic variable h, this process of itiner- 
ancy over different quasi-stationary states with dif- 
ferent spot configurations is understood as follows: 

Each quasi-stationary state corresponds to a state 
with a different configuration or a different number of 
spots, sp n - This state is also represented by a "sad- 
dle" point in the high dimensional dynamical systems 
of the collective coordinates h. There is a stable man- 
ifold to this saddle connecting from a state with a 
different number of spots, while the orbit spirals out 



of the saddle, corresponding to the breathing motion 
of a spot. Although this low-dimensional dynami- 
cal description of quasi-stationary state seems to be 
rather effective, the transition, indeed, occurs in a 
high-dimensional dynamical system. Such switching 
among effectively low-dimensional states within high- 
dimensional dynamical system is studied as chaotic 
itinerancy [341 |3"51 |4"!S] . The present model gives a 
novel example of chaotic itinerancy, whose mecha- 
nism has to be elucidated in future. 

When the number of spots is larger than or equal 
to 6, a variety of transition processes appear due to 
drastic increase in possible configurations of spots. 
As the spot number increases beyond 6, the com- 
binatorial explosion of the cofiguration of spots sets 
in, so that stable manifolds connecting many saddle 
points come to be entangled. This leads to a variety 
of transitions over quasi-stable spot states, resulting 
in chaotic itinerancy. 

In [131 HI] it is shown that combinatorial explosion 
which appears beyond the degrees of freedom 5 ~ 
9 leads to complex dynamics behavior with chaotic 
itinerancy. Complex dynamics for the spot number 
beyond 6 may be discussed from this viewpoint. 
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Appendix Derivation of evolution equa- 
tions 

Here we derive 2 component RD equation in terms 
of current density i, charge density <?'(= q + s), from 

First, (ll|^(p| are written as 

Jk - jk + \--J^j{(l k + Sfe ) ~ fe-i + Sfe -!) 

+{q k + s k ) - (qk+i + s k +i)}} =0 (1) 

(2) 
(3) 



d<7fc ds k 
jdik . qk n 



N 



I = Jk = const. 



(4) 



fe=i 



V = + + rJ k = const. (5) 

Now, from ([5]), we get 

N N 

NV = r Jk + r + Sfe ) = consL 



k=l 



k=l 



1G 



By transformation of variables q' k = qu + s k , and 



NV = rl 



c 



Q' = const. 



(9) 



is derived. 

Now we eliminate Jk,jk with ([TJ , ^ , Q to derive 
evolutional equation in terms of qk,ik- From ([l]),([2]), 
we obtain 

ik = Jk + -jyj{q'k-i + q'k+i - 2( 7fe) 



jk = ik + 



dt 



Accordingly. 



dt = Jk- ik + -^Q^k-i + q'k+i ~ 2 ik) 



(10) 



is obtained. From Jfe = -^r — ^» and v k = v(ik) in 
Fig|3ja), ([3]) and (10 1 are written as 



<K = V_q^ 
dt r rC 



dh _ q' k _ Sk_ _ v(ik) 
dt~lClC I ' 
while from Q, ^ and (10), we get 

N 



(11) 

(12) 



fc=i 



N 



1 

RC 



{q'k-i + q'k+i-H)} 



1 

RC 



(q'k-i + q'k+i -Zq'k)}- 



dt 

k=l 

N 

= E{** 

k=l 

When the boundary condition is Neumann or period 
one, / = X^fcLi *fc — const, is satisfied. Therefore, 
from ( 12 1, 



N ,. iV 



k=l 



Sfe w(ifc) 



z 



TV 



JV 



=>■ Q' = Sk + C u(ifc) = const. (13) 
fc=i fc=i 

is obtained. Since we have defined Sk so that it van- 
ishes when the value of each variables is spatially ho- 



mogeneous, (13 1 is rewritten as 



' = Cj2v(^) = Cv(jj)N = const 



k=l 



Hence Q',V is determined by ^, 



Q' = Cv( — )N = const. 

1 1 Q' 1 , 1 N 

V = r h — — = r hf — = const. 

N C JV N N 



got Q' = CE^i «(*'*) 



Here we discuss Sk that is defined as spatial charge 
distribution. According to Nasuno's experiment, 
macroscopic current- voltage characteristic V = V(I) 
is monotonous, therefore I and V are kept constant, 
whereas local current and charge density between 
electrodes are not so but can show spatiotemporal 
patterns. Hence we assume that local character- 
istic Vk = v(ik) has negative slope, as is typical 
in glow dischargee. Now, if there were not s k , we 
' ' '' ..(.-. 1 — const. In homogeneous 

case, (111, (12 1 actually satisfied 

CEk=lv(ik)(= Cfct, «(£) = Cq' k N) with- 
out Sfe. But, in general, ( fllj i, ( fl2| ) can't satisfy it 
without Sfc in inhomogeneous (i k ,q'k) case, and then 
need another term with physical meaning to satisfy 
it. Hence, we assume that some spatial charge dis- 
tributes besides q k - This is why we introduce Sfe. 

If Sk were independent variables, there would exist 
5JV valuables for AN equations. To solve the above 
equations, we need another constraint for the descrip- 
tion of Sk as a function F of {i m , <?m}, {jn =1,2, JV) 

Sfe = F({i m ,q m }) (m= 1,2,...,JV) 

This means that we assume each Sk reaches steady 
state quickly, that is, s k is eliminated adiabatically. 
This assumption appears valid, because Sk is charge 
density flowing in the longitudinal direction, and 
therefore the time scale of it is considered to be much 
faster than that of charge modulation in the perpen- 
dicular direction so that the former is separated from 
the latter, in the case of low resistive limit of elec- 
trodes (small but non-zero r). 

Hence Sk are determined so as to satisfy 

dik_ _ <ik_ _ £fe _ vjik) 
dt ~ IC IC I 

N 

I = ik = const. 
fe=i 

N I 
Q' = q'k = Cv(-)N = const. 

fe=i 



V = r jj + v (jj) = const - 



Now Sk is written as 



Sfe 



N 

|(Q'-C5>(4)). 



fc=i 



As understood easily, with this form of s k , it vanishes 
at the initial homogeneous state. By defining Sk in 
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this way, our equations are solvable. Considering spa- [21 
tial continuum limit and pressure effect as diffusion 
of i, the equations are described as 



dq' 


V 


q' 


~dt ' 


r 


rC 


di 


q' 


s 


dt 


IC 


~ Tc 



RC 



[22 

(14a) [23 



I 



DiAi 

s = j{Q' -C J v(i)dS} 
1=1 idS = const. 



[24 



(14b) 

^ [25 

(14d) I 26 



q'dS = Cv( — )S = const, 
o 

V = r— + v(-) = const. 



References 

S.Nasuno,Chaos 13,1010(2003). 



(14e) 



[27 
[28 

(14f) [29 
[30 
[31 



[10 
[11 

[12 
[13 

[14 

[is: 

[16. 

[i7: 
[i«: 

[19 
[20 



S.Nasuno, Report for Grant-in-Aid for Scientific Re- [32 
search(2001), (11837015), Kyushu Institute of Technol- 
ogy, unpublished. 

[33 

T.Sugawara,K.Kaneko,Prog.Theor.Phys.Suppl.l61,344(2006) 



V.S.Zykov,A.S.Mikhailov,S.C.Miiller, 
Phys.Rev.Lett.78,3398(1997). 

S.Grill,V.S.Zykov,S.C.Muller,Phys.Rev.Lett. 
75,3368(1995). 

G.Nicolis,I.Prigogine, 

Self-Organization in Non- Equilibrium Systems 
(Wiley, New York,1977) . 

T.Ohta,M.Mimura,R.Kobayashi,PhysicaD 34,115(1989). 

T.Ohta,Physica D 151,61(2001). 

S.-I.Ei,M.Mimura,M.Nagayama,PhysicaD 
165,176(2002). 

Y.Nishiura,D.Ueyama, Physica D 150,137(2001). 

Y . Nishiur a, T . Teramoto , K . I- 
.Ueda,Chaosl5,047509(2005). 



[34 
[35 

[36 
[37 

[38 

[39 
[40 

[41 



M.Suzuki,T.Ohta,M.Mimura,Phys.Rev.E52(4),3645(1995). I 42 



S.Kawaguchi,M.Mimura,SIAM 
J.Appl.Math.59,920(1999). 

A.M.Turing.Philos.Trans.R.Soc.London B,237,37(1952). 

P.Gray & S.K.Scott,Chem.Eng.Sci.38,29(1983). 

P.Gray &c S.K.Scott,Chem.Eng.Sci.39, 1087(1984). 

P.Gray &; S.K.Scott,J.Phys.Chem.89,22(1985). 

J.E.Pearson,Science 261,189(1993). 

K.G.Miiller,Phys.Rev.A 37,4836(1988). 

C.Radehaus,R.Dohmen,H.Willebrand,F.- 
J.Niedernostheide,Phys.Rev.A 42,7426(1990). 



[43 
[44' 



C. Radehaus,H.Willebrand,R.Dohmcn,F.- 
J.Nicdcrnosthcide,G.Bengel,H.-G.Purwins,Phys.Rev.A 
45,2546(1992). 

J.Berkemeier,T.Dirksmeyer,G.Klempt,H.- 
G.Purwins,Z.Phys.B 65,255(1986). 

Ch.Radehaus,T.Dirksmeyer,H.Willebrand,H.- 
G.Purwins,Phys.Lett.A 125,92(1987). 

Ch.Radehaus,K.Kardell, H.Baumann, D.Jager, and H.- 
G.Purwins, Z.Phys.B 65,515(1987). 

Yu.A.Astrov,Yu.A.Logvin,Phys.Rev.Lett. 79,2983(1997). 
Yu.A.Astrov.Phys.Rev.E 67,035203(R)(2003). 
Yu.A.Astrov,H.-G.Purwins,Phys.Lett.A 283,349(2001). 
Y. P. Raizer , GasDischargePhysics ( Springer, Berlin, 1991). 

D. D.Sijacic&U.Ebert,Phys.Rev.E 66,066410(2002). 

Y.P.Raizer,U.Ebert,&D.D.Sijacic,Phys.Rev.E70,017401(2004). 

Actually, this logarithmic characteristic is not essential 
but cubic characteristic only is enough for reproducing 
the phenomena. 

E. L.Gurevich, A.W.Liehr, S.Amiranashvili, and H.- 
G.Purwins, Phys.Rev.E 69, 036211(2004). 

In normalizing equations | |6a| ~|6f[l, we determined 
(i c ,v c ) so that the critical point equals to (1,1). 

K.Kaneko.Physica D 41,137(1990). 

K.Kaneko, I. Tuda, Complex Systems:Chaos and Be- 
yond, Springer(2000). 

L.P.Shil'nikov, Sov.Math.Dokl. 6, 163(1965). 

Here we display only the states generated by splitting, 
while complicated process involving collapse of spots is 
not included. 

E.L.Gurevich,Sh.Amiranashvili,H.-G.Purwins, 
J.Phys.D:Appl.Phys. 38, 1029(2005). 

Y.Kobayashi,M.Sano,Phys.Rev.E 73,015104(2006). 

K.Krischer,A.Mikhailov,Phys.Rev.Lett. 73,3165(1994). 
D.Battogtokh,M.Hildebrand,K.Krischcr,A.Mikhailov, 
Phys. Rep. 288,435(1997). 

N.B.Ouchi,K.Kaneko,ChaoslO,359(2000). 

K.Kaneko,Phys.Rev.Lett.78, 2736(1997). 
K.Kaneko,Physica D 124,322(1998). 

K.Kaneko,Phys.Rev.E66,055201(2002). 

S.Ishihara,K.Kaneko,Phys.Rev.Lett.94,058102(2005). 



18 



